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1. Introduction 



Abstract. A class of spherical collapsing exact solutions with electromagnetic 
charge is derived. This class of solutions - in general anisotropic - contains however 
as a particular case the charged dust model already known in literature. Under some 
regularity assumptions that in the uncharged case give rise to naked singularities, it 
is shown that the process of shell focusing singularities avoidance - already known 
for the dust collapse [ 10 1 - also takes place here, determing shell crossing effects or 
a completely regular solution. 
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The study of charged solutions to Einstein's equations is certainly an interesting 
topic since the beginning of Relativity, and the possibility for an astrophysical object 
to possess an appreciable amount of electric charge has been debated for a long time 

El H1I21- Anisotropic charged solutions were studied in 0, where relevance of charge 

5h h— h 

tkjrj effects is highlighted, and, under self-similarity assumption, in [1]. Many numerical 

studies on the subject exist - see e.g. (6) - even outside the realm of relativistic fluids, 

such as charged scalar fields IfToTl . 

In particular, charged collapsing object have attracted researchers' attention, for 

different reasons. On one side, if charge addition prevents the collapsing solution 

to become singular, Reissner-Nordstrom regular interiors could be modeled conse- 

quently (51 El- O n the other side, if a singularity forms, it is extremely intriguing to 

investigate the charge effects to the causal behavior of the solution, in particular if the 

singularity is completely hidden inside a trapped region or a naked singularity exists, 

in the latter case providing counterexamples to Cosmic Censorship Conjecture. 

In ifTOll the charged spherical dust cloud collapse was studied to find that shell fo- 
es 

cusing singularities are avoided in general, obtaining either shell crossing singularities 
or completely regular solutions, depending on mass-to-charge ratio. In the same pa- 
per it was conjectured that shell crossing effects could be averted if spatial gradient of 
pressure remains nonzero. In the present paper we determine analytically a class of 
charged spherical exact solutions undergoing gravitational collapse, and study their 
final state. These solutions are in general anisotropic, but contain the charged dust 
model studied in IfTOll as a limit case. We find that shell focusing singularity forma- 
tion is avoided for non central shells even for this wider class of solutions, and the 

solution either remains completely regular, or shell crossing singularity is developed. 
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The paper is organized as follows. Section [2] is devoted to introducing the model 



in particular, in subsection 2.1 the procedure of charge addition to a spherically 



symmetric solution is sketched, while the class of (uncharged) background models 



is introduced in subsection 2.2 - see also the Appendix for the full list of conditions 
satisfied by these models. The new (charged) solutions are derived in Section [3] and 
their relevant properties are studied in Section |4j whereas conclusions are left for the 
final section. 

2. Charged spherical elements 

In this section we introduce relativistic charged models in general, and recall the 
class of model that we will perturb in next section using an electromagnetic charge. 

2.1. Addition of an electromagnetic charge. We start from a spherical material de- 
scribed by the line element, written in comoving coordinates (t,r,9,(f>), 

(1) ds 2 = -e 21 ^ dt 2 + e 2X ^ dr 2 + R(t, rf dfi 2 , 

carrying a non-vanishing electric charge density p e . We introduce the skew-symmetrc 
Maxwell tensor F^ v and, due to spherical symmetry and the absence of magnetic 
monopole iTTOl [T2| . the only nonzero component is F 01 . From the Maxwell's equa- 
tions 



(2) F"*.„ = 47rp e u", F^ v . p] = 0, 
being u a timelike unit vector field, we find that 

(3) F 01 = ^ 

where Q = Q(r) is arbitrary and is related to the charge density by the relation 

(4) Q'(r) = A7T Pe e x R 2 . 

Let us consider Einstein equations G = 8nT for the line element ([T]), where T = 
^mat + T em is the sum of the contributions of the stress energy tensor induced by the 
material structure and the stress energy electromagnetic tensor 

(5) (T aP ) em = ~(F£Ftf + ^g a pF^Fn, 
that using ([3]) is found to satisfy 

and therefore is diagonal when written as a (1, l)-tensor. The above procedure allows 
to add a charge to a given spherical model. 
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2.2. Spherical elastic-solid models. In the following we are going to describe the 
spherical model that will be perturbed with an electromagnetic charge following the 
above scheme. The model considered here has been found and fully described in 
B71[8]|, but for sake of completeness we briefly recall here the basic properties. 

We consider a source of gravitational field given by an elastic material under isother- 
mal conditions. The stress energy tensor in comoving coordinates is diagonal, 

(V) (T^ Q ) mat = diag (-e(t, r),p r {t, r),p t {t, r),p t {t, r)) , 

where e is the energy density and p r and p t are the radial and the tangential pressure, 
respectively. As described in ifTTl 1711. the additional information to close the system 
of Einstein equations are supplied by a constitutive function w = w(r,R,r)) where 
r) = e~ x . Note that w is a function of the spatial coordinates and of the space-space 
part of the metric, not depending on the angular coordinates 9 and for symmetry 
reasons. The quantities in (|7|) are related to w by the identities 

dw 1 dw 

(8) e = pw, p r = 2pr] — , Pt = -- p R— i 

where 

Vv 



p 



8irE(r)R 2 

is the matter density and E(r) is an arbitrary function. 

This framework suggests, as a more convenient way to describe the metric, the use 
of area-radius coordinates (r, R, 9, 0), in such a way that the line element becomes 

ma 

(9) ds 2 = -A(r, R) dr 2 - 2B{r, R) dR dr - C(r, R) dR 2 + R 2 dfi 2 . 
In [El [8]| it is described how, if we choose a r] _1 / 2 -linear constitutive function, 

(10) w(r, R, rj) = E{r) + *^(r, R)v~ 1/2 

determined by the choice of two free functions \l/(r, R) and Y(r, R) - constrained at 
most to fulfil some reasonability conditions that we will review in a moment - then 
Einstein field equations for the metric d9j) takes the form 



(ID ds 2 = - ( 1 - ^) G 2 dr 2 + 2G -dRdr - -dR 2 + R 2 dtf, 

R u u 



where 



(12) M 2 =y z += __ 1) 



2^ 
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and the function G is given in terms of a quadrature: 
(13) G(r, R) = f (-) (r, a) da + 



R Y(r,a) \u J r ' Y(r,r)u(r,r) 

It is quite straightforward to verify that the arbitrary function \P is the Misner-Sharp 
mass = f (1 — fl^V-R, Vi?)), whereas the Y is related to the comoving description 
of the matter through the relation Y = R'r). Also observe that uGr] = 1 as one can 
find translating the line element ([TJ into ( [TTj ). 



Using and ( fT0[ ), the components of T mat as functions of (r, R, rf) are given by 



(14) e = - n9 ' r , - + 



4irR 2 YuG AttR 2 ' 

1 j Pr ~ 47T.R 2 ' P * ~ Stt^G \Y J r 8ttR 

The above class contains, as particular cases, Tolman-Bondi-Lemaitre dust models - 
corresponding to ^ = \&(r) and Y = Y(r) - and vanishing radial pressure models - 
taking ^ = \I/(r) and letting Y possibly depending on R too. In [8 J some particular 
cases separately studied in literature including the above cited are reviewed. 

As said before, in order to describe reasonably a collapsing matter, some conditions 
must be placed on the arbitrary functions \&(r, R) and Y(r, R). Using the notation 
from [7, Definition 2], we suppose that they are chosen in such a way that ( [TTj ) is a 
collapsing area-radius separable (ARS) spacetime - see Appendix |A| for a complete 
list of assumptions. This implies, among the others, that w.e.c. is satisfied and that \& 
and Y are positive and regular functions whose Taylor development in a neighborhood 
of (0, 0) is given by 

(16) *(r, R)= Yl ^ riRJ + 2 + ° 3+m(r ' ^' 

j+.j=3 i+j=3+m 

(17) Y(r, R) = l+^2 ViJ riRj + ^ riRj + ° 2 +^ r > 

i+j=2 i+j=2+£ 

Note also that in view of the Definition given in Q it must be ^30 > and ^ r , ty R > 
0. 

Finally, in order to obtain a global model, a matching with an external space will 
be performed at £ = {r = r^}, requiring that the fundamental forms of the two 



metrics at £ coincide (Israel-Darmois junction conditions). From ( |T5j ) we observe 
that the radial pressure p r in general does not vanish along S, which as well known 
is a necessary and sufficient condition to match the solution with a Schwarzschild 
exterior. In this ore general case, a natural choice for the exterior metric is given by 
the generalized Vaidya spacetime [[TBI 

(18) Sext = - (l - dV 2 + 2dVdS + SW, 
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where fi(V, S) is an arbitrary (non negative) function. Reissner-Nordstrom metric is 
contained here as a particular case (p = p — with p^ constants). 

The immersion of S in the two spacetimes can be parameterized respectively by 

(a,6,4>) (r&, cr, 6,(f)) and (a,9,(f)) ^ (V(a), S(a),9,(j)), and junction conditions 
are found to be 



(19) S(a) = a, 



dV(a) 1 



dcr u(Y — u) 



\{r b ,R=a) 



(20) fi(V(a), S(a)) = tffo, a), J£(V(*), S(a)) = 0. 

3. Charged elastic models 
Let us consider the spherical model above recalled and let us add a charge density 



as described in Section 2.1 Coupling between elasticity and charge gives rise to a 
field theory where Lagrangian density of the elastic material ([TT| A mat = —yj—ge = 
— y/^gpw - where p denotes the matter density, p = rj/ 1 (AiiE{r)R 2 ) - sums up 
with the Maxwell Lagrangian density A em = —\^/—gF^F^. Then using Belinfante- 
Rosenfeld theorem the total stress energy tensor T is now given by the sum T mat + T em , 
that in view of ([6]) and |7]) is diagonal when expressed in comoving coordinates. More- 
over, the perturbation given by ([6]) depends on (r, R) only. This suggests the idea that 



the new model may belong to the class considered in Section 2.2 In particular, if we 



search for new functions ^i, Y\ such that the new stress energy tensor obeys rela- 



tions ( [14] ) and ( [15] ) for these new functions, it is easy to check that we are led to the 
following definition: 

Q{r? 



(21) ®i(r,R) = *(r,i?) - 



2R 



In other words, the information about the charge effect can be completely encoded 
into the mass and the Y function in such a way that the solution obtained represent 
a fictitious material belonging to the same class, though with a new choice of the 
arbitrary functions. Note that this procedure is possibile since charge in ([6]) is added 



to the off-shell relations ( fl4] ) and ( [13] ) where 77 is not explicitly depending on (r, R). 

Of course, the new degree of freedom introduced with the charge modifies the con- 
stitutive relation between the elements of the total stress energy tensor, and indeed 
now we have 

(23) Wl (r, R, v) = w(r, R, 77) + ^r?T 1/2 

where w(r, R,r/) is given by ( [TO] ). At this stage, however, a remark is in place: al- 



though \P and Y satisfy all the properties sketched in Section 2.2 and listed in the 



Appendix - especially regularity (16) and ( 17 ) -in view of (21 ) and (22) one cannot 
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expect that the same properties are satisfied by \&i and Y\. Put it in another way, 



choosing the two arbitrary functions as in ( |2T| ) and ( |22| ) above, the charged model ob- 
tained is not an ARS collapsing spacetime in the sense of [|7l Definition 2], and then 
results from there do not work here. For instance, we cannot exclude shell crossing 
singularity formation - that on the contrary is likely to occur as we will see later. Nev- 
ertheless it is easy to see that the new model automatically satisfies the weak energy 
condition if the old does. Moreover, if we require regularity of the energy density in 
the center, recalling that R = r at initial time, the condition Q 2 (r) = 0(r 4 ) as r — > 



is needed, as it can be seen using ( fT4| ) together with the facts that R' = uGY and 
R' = 1 at initial time. Assuming regularity of the charge function Q(r) as done for 
the other arbitrary functions, this amounts to say that Q(r) = q r 2+p + o(r 2+p ), with 

P > o. 

Moreover, we have to assume that local isotropy at the center continues to hold 
once the charge is added. In comoving coordinates, using ([6]), that would mean that 

(24) Urn (p rl (t,r) -p n (t,r)) = lim [p r (t,r) ~Pt(t,r) - = 0. 
Supposing that the local isotropy holds for the uncharged model, this means that 

r ^o+ R 4 (t,r) 

whenever we tend to the regular center. Recalling again the initial data R = r, that 
implies that Q(r) = o(r 2 ), and therefore the condition stated above on the behavior 
of the charge function must be strengthened as follows 

(25) Q(r) = q r 2+p + o(r 2+p ), p > 1. 

4. Collapsing solutions and singularity formation 



In view of the lack of regularity at the center - see plj ) and ( |22] ) - one cannot expect 
that these solutions collapse completely to a shell focusing singularity. To begin, an 
important preliminary result can be immediately stated. 

Proposition 4.1. Shell focusing singularity cannot form from the collapse ofnon cen- 
tral shells of the charged elastic models described in Section^ unless a shell crossing 
singularity forms at an earlier comoving time. 

Proof. We will actually show that, for each r > 0, the function R !->■ Yi(r , R) must 
vanish for some R e]0, r [. Since the equation Y 1 = controls shell crossing singu- 
larity formation, we will deduce that shell focusing non central singularities cannot 
form before possibile shell crossing singularities. 

From ( |22| ), and recalling from ( [T7| ) that Y(r,R) is strictly positive, the zeroes of 



Yi corresponds to zeroes of ^x, r (r, R) = \l/, r (r, R) — Q(r ^ ^ , Fixed r = r , this 
quantity negatively diverges as R — > + , which means that the continuous function 
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R i — y \l>i, r (r ,-R), that starts positive for R = r , must vanish for some positive 
R<r . □ 

In view of the above result, non central shells either develop a shell crossing singu- 
larity or they stay regular, collapsing indefinitely to a stable non singular state or even 
bouncing and finally entering a period of expansion. 

Shell crossing formation corresponds to the zeroes of the function Y (that in comov- 
ing coordinates reads as R'rj), while no singularity formation is possible if ^ — l+Y 2 
vanishes - consider ( fT2[ ) restated using ^ and Yi in place of Y . 

Actually, points of the (r, R) plane where ^ — 1 + Y x 2 is negative describe a region 
where the dynamics of the spacetime is not allowed; taking initial data such that this 
quantity is positive, and studying the evolution of this quantity for a fixed shell r = r , 
if the quantity vanishes for some R = R u (r ) and attains negative values for smaller 
values of R, this means that the shell labeled r cannot collapse beyond the limit 
value R = R u (r ), then indefinitely tends to it or arrives there in a finite amount of 
comoving time, bouncing back and entering an expanding phase. 

In view of the aforesaid, there exists a curve R sc (r) such that, for r > 0, R sc (r ) 
is the greatest zero of the function R !->■ Y%(ro,R) between and r (recall that 
Yi(r , r o) > by the initial data choice). Recalling ( |2~2~| ) for Y\, we have to determine 
the zeroes of (^i), r and using ( [To] ) and ( f25] >, it is found explicitly that 



(26) R sc {r) = ( 2 +^°> +i + o(r ^+i). 

3^30 

Moreover, to investigate the physics of the collapse, we need to study the apparent 
horizon, implicitly defined by solution of equation 

(27) R = 2^ 1 (r,R). 



Recalling equations ( |16| ) and ( [25] ) we find the existence of two different horizon curves 
satisfying ( |27| ), whose behavior depends on the charge function Q(r). In particular, 
we have the following situation: 

• when p > 1 (and ip 30 > %)'■ 

(28) R hl (r) = 3^r 2 v +l + o(r> +1 ), R h2 (r) = 2^ 30 r 3 + o(r 3 ). 

^30 

• when p — 1: 

(29) 

Rhi(r) = Urn - yfifaf-qo 2 ) r 3 +o(r 3 ), R h2 (r) = U 30 + \j ^ - qA r 3 +o(r 3 ) 



As can be inferred from ( 29 ), when p = 1 and ip 30 < qo,no horizon forms. In all other 



cases, we have two apparent horizons which bound the shell crossing curve R sc (26) 



8 



F. CIPOLLETTA, R. GIAMBO 



one from below and the other from above. As said before, vanishing of the function 



(30) Ul (r, R) := y Wli ^ R) - 1 + Y 1 (r, Rf 

could prevent shell erasing formation. From the above analysis, there is certainly a 
case when this happens, that is when no horizon formation takes place. Indeed, in that 
case it must be R > 2^i throughout the whole evolution, and then on the curve R sc 



( 26 ), the quantity under square root in ( 30 ) is negative. Therefore, u\ already vanished 



before R sc which means that the solution remains regular for all times. To gain more 
information on the behavior of u\ in this case (p = 1, -?/>30 < <?o)> one can study the 
quantity Ru\ along curves (r, kr) for k e]0, 1], finding that 

R Ul 2 (r,R) lR ^ kr = V (ky + o(r 3 ), 

where 

2A(k) 

(31) V (k) 



(3^0 + 2fc^ 21 + k^ 12 ) ' 
with 

(32) A{k) = -3q 2 + (A^i2 + 2^ 2 i + 3^ 30 ) • 

• (k 3 ip Q3 + k 3 y 02 + k 2 ipi2 + k 2 y u + kip 2 i + ky 20 + ^30) 



The denominator in ( pTj ) is positive - is the leading term of \& ir (r, kr) which is positive 



by definition, see subsection 2.2 - therefore we have to study the sign of A(k) < |32] ) . In 
particular we have A(l) > (as it is the initial data) and A(0) < (since ^30 < Qo)- 
Thus (f(k) changes sign for some k e]0, 1[, and then the zeroes of Ru\ defines a 
curve 

R u {r) = k r + o(r) 

which bounds R sc from above as expected. The collapse behavior is sketched in 
Figure [TJ where a graphic is plotted with R/r on the vertical axis to get rid of the 
central coordinate singularity. The collapse starts from R/r = 1 and every shell 
collapses until it reaches R = R u (r) where it cannot proceed further since the quantity 
under square root in ([30]) would be negative. Then, once the shell radius reaches 
the value R = R u (r), it starts bouncing back entering an expanding phase. Since 
Ru( r )/ r — > ko > as r — > + , even the center remains regular for all times. 

Remark 4.2. It can be seen that the quantity A{k) is given by the leading order term 
of the function ^, r {r, R)Ru 2 (r, R) - 2Q(r)Q'(R), evaluated for R = kr. 

Now let us examine the other cases. When p = 1 and ^30 > Qo we can perform a 
similar study as before, finding that Ruf(r, kr) = Lp(k)r 3 + o(r 3 ), with (p(k) given by 



pTj ) above, and then its sign is related to the sign of A(k) ( |32| ) again. However, in this 
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FIGURE 1. When a curve R u (r) - where the quantity under square 



root in pO) changes sign - is defined, there exists a region - the shaded 
area below that curve - that cannot be reached starting the collapse 
from R/r = 1. 



case both A(0) and A(l) are positive, so we have to consider separately the following 
two subcases: 

(a) if A[k) admits one root in the interval ]0, 1[, then as before it remains defined 
a curve 

Ru(r) = k r + o(r) 

which bounds R sc from above, and then the situation is the same as that de- 
picted in Figure [T] 

(b) if A(k) > in [0, 1], then we have to evaluate Ru\ along R = kr 2 . For sake 
of simplicity we consider here only the case ^30 > qo (the special case ^30 = 
go can be discussed with similar methods). It can be found Ru\(r, kr 2 ) = 
Lp(k)r 3 + o(r 3 ), where 

2(^so 2 - go 2 ) 



(33) 



(34) 



<p(k) 



.30 



that is positive. Evaluating finally Ru\ on curves R = kr 3 , one gets again a 
behavior of the kind Ruf(r, kr 2 ) = Lp(k)r 3 + o(r 3 ), where 



<p{k) 



2(^30 2 - Qo 2 ) 



+ 



qo' 



(qo 2 - ^3o 2 )- 



^30 ' ktjj 30 

and therefore the zeroes of Ru\ describe a curve R u (r) whose behavior is 
given by 



RJr) 



^ r 3 + (r 3 ). 



30 



Comparing the above equation with ( f26| ) we conclude that R u (r) < R sc (r 
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Figure 2. When the curve R u {r) lies below the inner horizon, then a 
shell crossing singularity - denoted by the curve R sc (r) - forms. The 
shaded area represents the trapped region. 

for r > and then in this case a shell crossing singularity develops just after 
the formation of the outer horizon. The situation is represented in Figure [2} 

Finally, case p > 1 is yet left to be studied. It is easy to check that zeroes of Ru\ 
corresponds to solutions to the equation 



2^ r 2 R u + R U 2 ^J [-1 + Y 2 } - 2QQ'R u Y 2 V yr = Q 2 



*, r 2 - Q' 2 Y 2 



From this, using expansions ( [16] ), ( fr7| ) and ( |25| ) and comparing the leading term of the 
left side member with the leading term of the right side member in the above equation 
leads to the following behavior for function R u such that «i(r, R u {r)) = 0: 

a, 2 

(35) R u (r) 



2if>. 



30 



An immediate comparison with equation ( [26] ) shows that R u (r) < R sc (r) for positive 
r and then all shells as before undergo a shell crossing formation, and the situation is 
again represented in Figure |2j 

Summarizing, we have that the behavior of the collapse is driven by the parameter 



p defined in (25 ) in the following way: 



Theorem 4.3. Under the assumptions made in Section^a shell crossing singularity 
forms up to the centre when: 



(1) the parameter p defined in ( [25] ) is greater than 1, or 

(2) p — 1, > go and the function A(k) defined in ((32]) is positive VA; G [0, 1]. 

In all other cases, provided ^30 7^ gcQ^ e solution remains regular, and no singularity 
forms up to the center. 



the highly non generical case p = 1, -030 = Qo can be treated with similar methods as before, just 
taking into account higher order terms of the arbitrary functions. 
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5. Conclusions 

The class of spherical models carrying an electromagnetic charge considered in 
this paper has been explicitly derived exploiting the description in terms of a state 



function w = w(r } R, 77) - see subsection 2.2 Since the addition of a charge modifies 



the stress-energy tensor with a quantity depending on (r, R), it is possible to derive 
the state function for the new solution, seen as a fictitious material. Therefore, one 
can say that the process of solution charging is - under a mathematical point of view - 
internal to the class of models considered, i.e. those determined by the state function 
( flQ] ). Although the charged solution does not satisfy all the regularity properties of the 
original one, however it is possible to explicitly derive the metric and then to study its 
properties. 

As we have seen, the electromagnetic charge deeply changes the final state of the 
solution. It is a well established fact indeed that the uncharged solutions here consid- 
ered always develop a shell focusing singularity that in some case is naked - depend- 
ing on Taylor expansion of the arbitrary functions \l/(r, R) and Y(r, R), see Em for 
full details. Instead, here the formation of a shell focusing singularity is impossible 
for non central shells, and otherwise either the solution bounces back and enters an 
expanding phase, or a shell crossing singularity takes place. As is well known shell 
crossing singularities are Tipler-gravitationally weak [fi3l and the spacetime may ad- 
mit an extension beyond this kind of singularity, possibly developing a shell focusing 
singularity at a later comoving time. This issue still remains an open problem in gen- 
eral, that has been partially solved only for particular cases, as Tolman-Bondi models 
[fi4ll . where metric can be extended in a natural but non unique and, yet bounded, non 
regular way beyond the shell-crossing singularity. 

The occurence of the two situations - either regularity or shell crossing effect - is 



determined (Theorem 4.3 ) by the Taylor development of Y and the charge function 
Q(r). Notice that, unlike the charged dust cloud case studied in ifTOl . the final state 
is not determined exclusively by the charge-to-mass ratio at initial time t = 0, since 
now the mass depends on R also and at initial time (when R = r) we have ^(r, r) = 
(^2i+j=3 V'yJ 7-3 + o(r 3 ), that in general does not allows to determine the coefficient 



^30 needed in Theorem 4.3 



However the avoidance of shell focusing singularities - up to spacetime extension 
as said before - already found in ifTOll remains true for this wider class of (generally 
anisotropic) solutions, giving a first answer (in the negative) to the conjecture raised in 
[fTOll that the occurence of shell crossing singularities was determined by the vanishing 
of the pressure spatial gradient. It emerges that electromagnetic charge completely 
changes the final state of the collapsing star in this anisitropic case. 
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Appendix A. Collapsing area-radius separable (ARS) spacetimes 



The models introduced in B [H and recalled in Section Z2] completely depend 
on the choice of two functions \& and Y. An ARS spacetime is described by the 



metric ( fTT| ), where (fT2])-(fT3~]) holds, and \I/(r, R), Y(r, R) are two regular functions. 
Since we consider initial data such that R = r, the spacetime is defined in the set 
S = {(r, R) : < r < r^, < R < r}, and defines the junction hypersurface 
with the external solution - see end of Section [2] above. Now we give a complete list 
of the properties that these functions must satisfy to obtain what in [0 is termed as a 
collapsing ARS spacetime. 

• the arbitrary functions must be positive up to possible singularity formation, 
*(r, R) > and Y{r, R) > 0, when R > 0; 

• the weak energy condition (w.e.c.) must hold. A sufficient condition is given 
by 

(36) M>, r > 0, * )R >0, *,r>f^(^J , *,r>**,r R ] 

• regularity and isotropy of the metric at the center of symmetry must be im- 
posed. It can be seen that this is equivalent to the following conditions: 

(37) F(0,0) = 1, 

(38) (0, 0) = , r (0, 0) = ^(0, 0) = * rr (0, 0) = * >rR {0, 0) = # ^(0, 0) = 0; 

• the initial density is required to be decreasing outwards, which means that 



(39) * rr (r, r) + 2^, rfl (r, r) + ^ jRR {r, r) - - r (r, r) + # R (r, r)) < 0; 

r 

• shell-crossing singularity formation must be excluded. This leads to the (suf- 
ficient) conditions 

(40) ®A r , R ) +Y(r,R)Y >r (r,R)R> 0, 



(41) 



o Y(r,a) \u 



(r, a) da 



u(r, r)Y(r, r) 



> for r > 0; 



• shell-focusing singularities must form in a finite amount of comoving time. 
The mathematical conditions to express this requirement are 

(42) K r (0, 0) = Y R (0, 0) = 0, (* >rrr (0, 0), * rrJi (0, 0), V >rRR (0, 0)) + (0, 0, 0). 



In particular, regularity of \& and Y, together with (]37[)-([38|) and ( [42] ) allows to write 
developments ( fT6] ) and p7| ). Moreover in [[7l [8]| to ensure the very last condition in 



( |4"2~1 ), it is assumed that ^ )r . rr (0, 0) > 0, that in view of ( fT6] ) means > 0. 
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